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Properties of X(3872) are studied by regarding it as a DD* hadronic molecule with J PC = 
2 _+ in the phenomenological Lagrangian approach. We find that our model with about 97.6% 
isospin zero component explains the existing data nicely, for example, the ratio S(X(3872) — > 
J/ipiT + iY~ it ) /B(X (3872) — » J/ipTv + iv~). We predict the partial widths of the radiative decays 
of X(3872) jJ/ip, jip(2S) and the strong decays of A(3872) ->■ J/ip-K + TY~, J / ipn + n~ it as 
well as X(3872) — > XcJ^°- Our analysis shows that the measurement of the ratio Z3(X(3872) — s> 
Xco7T°)/B(X(3872) -> Xcitt ) may signal the nature of X(3872). 

PACS numbers: 14.40.Rt, 13.20.Fc, 13.25.Ft, 12.60.Rc 



I. INTRODUCTION 



A" (3872) is one of the mysterious exotic states which was observed by the Belle Collaboration in the exclusive 
B ± -> ^Tr+Tr" J/i/j decays in 2003 Q. Later, this state was confirmed by the CDF jl, the DO |[ and the BaBar 0- 
0] Collaborations. The observation of its decay channels 'yJ/tp [E and t4>(2S) jl| indicates that A"(3872) has a 
positive charge parity. By analyzing the invariant 7r + 7r~ mass distribution Q and the total angular distribution [l(3] 
of the J/t/jir+ir^ decay mode, the CDF Collaboration has narrowed its possible spin-parity (J p ) to be 1 + or 2~. 

The X(3872) is identified as an "exotic" state because its observed mass does not fit into the quark model predictions 
for the charmonium states pdl4l3j and, the branching fraction of X — > Jipui — > J/ipn + ir~ir and X — > Jipp — > 
J/ipTr + n~ is almost equal. In the literature, A"(3872) has been interpreted as a ccg hybrid state [3 EH, a diquark- 
antidiquark state [l6| , a cc and meson cloudy mixing state [l7-20] or a deuteron-like molecular state concerning its 
mass is approximate to the threshold of D°D*° [2l| - (34| . In addition, this state has also been studied with the coupling 
channel approach [35l - l37j . For other studies of X(3872) with J PC = 1 ++ , see the recent review Refs. [H, and the 
references therein. 

Recently, the analysis of the X(3872) — > J/tpuj — > J/ipn + ir~Tr decay mode by the BaBar Collaboration 0] indicated 
that X(3872) favors a negative parity, i.e., J PC = 2 ^. Then, the radiative decays of X(3872) — > jJ/ip, r )ilj(2S) 
were studied by assuming X(3872) with J PC = 2 ^ as the r/ C 2(lD) and it was found that the data cannot be 
interpreted. This indicates that X(3872) may be an exotic state. 

In this paper, we study the decay properties of A" (3872) in the effective Lagrangian approach by regarding it as a 
DD* molecule with J PC = 2 '". We use the method which was proposed by one of us in the previous work [4li - [43| 
and was applied to study the properties of A"(3872) with quantum numbers J PC = 1 ++ [H, HI] and other exotic 
states pil - lia ]. The composite state X(3872) is defined via the compositeness condition Z — which was originally 
used in the study of deuteron [4?], |48| which is a bound state of proton and neutron with Z = as the wave function 
renormalization constant of the composite particle. 

We write the wave function of X(3872) as 

|X(3872)) = C -^\D°D*°) + ^\D + D*-) + Cc, (1) 

which is in terms of charge eigenstate. Or equivalently, one may regard the observed X(3872) as a mixing state of 
1 = and 1=1 states, i.e., 

|X(3872)) = cos <£|X(3872)) /=0 + sin </>|X(3872)) J=1 , (2) 
with cos 9 = (cos <f) + sin 0) /y/2 and sin 9 = (cos (f> — sin <fi) / \[2 and 

|X(3872)) J=0 = l(\D°D*°) + \D+D*-)^j + Cc, |X(3872))/=i = ^(\D a D*°) - \D + D*-)^ + Cc. (3) 

We find, using the two parameters, the mixing angle <fi between the isospin singlet and triplet components introduced 
in Eq. (|2|) and the size parameter Ax introduced in the following to illustrate the distribution of the constituents 
in the molecule, we can explain the date quite well. From the magnitude of the mixing angle <p which was fitted 
from the data, we see that, in the wave function of A" (3872), the isospin singlet component is dominant which is 
consistent with the analysis from the phase space. That is, in case of X(3872) as a tensor meson, the phase space 



of the decay A(3872) —> J / '^jm + n~ n° is about one percent of that of the decay A(3872) — > J/ipn + Tr~ so only few 
percent of isospin one component in A (3872) can accommodate the large isospin violating partial width. We also 
explicitly calculated other decay widths of A(3872) — > jJ/ip, r yi(j(2S) and A(3872) — > XcjTi" («/ = 0, 1, 2) using these 
two fixed parameters. Comparing our results with the more precise data observed in the future, one can determine 
the structure of A (3872). For example, if the future data are not consistent with our present results, there may be 
J/ipto, J/ipp or/and cc components in the A(3872) wave function. 

This paper is organized as follows: In Sec. [ill we propose the theoretical framework applied in our calculation. In 
Sec. IIIII we give the formalism for the strong and radiative decays of A(3872) explicitly. Sec. IIVI is our numerical 
results and discussions and the last section is our conclusion. 



II. THE THEORETICAL APPROACH 



As was mentioned above, we regard A(3872) as a DD* bound state with J PC = 2 + with the explicit wave function 
((T|) . Following Ref . |49| , one has the free Lagrangian for tensor field as 
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^ V" T~\fJtV]\<J y 

free — « ^ -^-Aff) 



(4) 



where the symmetric tensor X^v = X ufl denotes the J PC = 2 + field for A (3872) and 

jyiv-M = (n + ml){-(g l *g'" r + g" x g^) - g^g Xa ] 

+g x<T d f *d v + g^ v d x d' y - ~ (g va d^d x + g vX d^d^ +(fj,^u) 

with m x being the mass for A" (3872). The propagator for A AI „(3872) is obtained as 

d^p i 
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Taking into account the quantum numbers of A (3872) and its constituents DD*, one can write down the effective 
Lagrangian describing the interaction between A (3872) and its constituents as 



£ x = -^=X liV (x) J dxidx 2 ^x{{x\ - x 2 ) 2 )8(x - uj v x\ - lo p x 2 ) 

1 



(7) 



x {.<?* [cjLixu x 2 ) + C^(x u x 2 ) - -g^C^ a (x u x 2 )] + g c x [c°,{x u x 2 ) + C^{x u x 2 ) - -g flv C^ a (x 1 ,x 2 )] }, 

where g x is the effective coupling constant for the interaction between A(3872) and its neutral constituents (D°D*° + 
D°D*°) while g x is for the interaction between A(3872) and its charged constituents (D + D*~ + D~D* + ). uj v and 
oj p are mass ratios with definitions 



uj v = 



(8) 



where m D and m D , are masses of the constituents pseudoscalar meson D and vector meson D* , respectively. In the 
Lagrangian (J7J, to illustrate the finite size of the molecule, the function &x((xi — x 2 ) 2 ) with the Fourier transform 



®x(y 2 



d 4 p - 



(27T) 



(9) 



has been introduced. And, for simplicity, we have defined the tensor C Mt , as a function of the constituents with the 
explicit form in terms of the neutral constituents as 



C» v {x u x 2 ) = D^{x l )d u D\x 2 ) + Dl\x l )d^{x 2 ). 
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Substituting the fields D° and D*° with the corresponding charged ones, one can get the explicit form of C^ v . 

Relation between the mixing angle 9 defined in Eq. ([]} and the coupling constant g x (g x ) can be yielded with the 
help of the compositeness condition Zx = with Zx as the wave function renormalization constant of X(3872). Zx 
is defined as the residual of X(3872) propagator, i.e., 



where g x H x (p 2 ) relates to the mass operator via the relation 

nf ;Q V) = J(r^ + /V a ).9^ x (/) + -, (n) 

with " • • • " denoting terms do not contribute to the mass renormalization of X(3872). The mass operator of X(3872) 
is illustrated in Fig. [T] Explicitly, we have the following relations between the mixing angle 9 and g N and g c 




FIG. 1: The mass operator for X(3872). 



cos 9 



(12) 



where T, x corresponds to the case with neutral constituents but corresponds to the case with charged constituents. 

The mass operator can be calculated from the Feynman diagram depicted in Fig. Q] In our calculation we applied 
the Gaussian form of <3>x(p 2 ), i.e., 



*jc(p 2 ) = exp(p 2 /A2), 



(13) 



with the size parameter A x introduced to parameterize the distribution of the constituents inside the molecule. In 
the following calculation, we will take A x as a free parameter and fix it from the data. 

Because of the uncertainty of the X(3872) mass measurement, we express the mass of X(3872) in terms of the 
binding energy as 



(14) 



where AE > is the binding energy and m D ,o = 2006.97 MeV and m D o = 1864.84 MeV [50]. In the following 
calculation, we take AE — 0.5, 1.0 and 1.5 MeV. 

To calculate the strong decays, we will use the effective Lagrangian including the following terms 
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+l g D , D , p {D;(D:d»T ■ n + d;(t ■ pJ»D*») + T-p»{Dld»D* v )l 
-i9DDtp n ipn; n(Dd M D) - ignn^oj^Dd^D) - ig D D P Pp. • (Drd^D), 

+g D *D^ vafi u^D* a pD + g D 'D P e^ p^ ■ D* a0 rD + H.c, 
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£>XcoDD = 9 X coDDXcoDD, 

C Xe0 D*D* = 9 X caD*D*XcoD*D* IJ -, 

C XclD *D = ig X ciD*DXcl;fiDD*^ + H.C., 

ZxczD-D* = 9x^D*D*Xc2^v[D*»D* V + D* v D*% (15) 

where ip n denotes J/ip and ip(2S) vector mesons, 7? and p stand for the pion and rho meson triplets respectively and 

D = (D a ,D + ) T and D* = (D*° , D* + ) T are the pseudoscalar and the vector charmed meson doublets respectively. 

♦+ ++ 

V^v (V — ip n , D* ,ui, p) is defined as V^„ = d^Vv — d v V^. The definition of the right-left derivative Ad B is Ad B = 
A(dB) — (dA)B. In Ref. [5l|, the coupling constants gD*D*v an d goov were determined with the help of the VMD, 
i.e., 

gvD'p = gDD P = -2.52, gD*D*Lu — gDDu — 2.84 , gD*D*^ — gDD*p — 7.64 . (16) 
Using the VMD and the partial widths T(J/<ip -> e+e") = 5.55 KeV and T{ip(2S) -> e+e~) = 2.38 KeV H3, we have 



gn^* ^g_oo^^ = L82 (17) 

gD*D*4> gDDip m$ y 1 (i/> — > e+e ) 

which leads to gw d*^' = goo^ 1 — 13.90. 

Prom the HHChPT including vector mesons within the hidden local symmetry method (52T45^ | (see Ref. [HH for an 
equivalent approach), one can yield the following relation for the coupling constants [54 [ 



A / m D , 

gD'Dui = gD*D P = 77ffW ' (1°) 

2 V m D , 

where g is the universal coupling constant introduced in the hidden local symmetry method with g = 5.8 ±0.91. And, 
A gives the coupling of the light vector meson with the heavy states and as analyzed in Ref. [HI], A = —0.41 GeV -1 . 
Finally, we get g D 'Du = 9d*Dq = -1-23 GeV -1 . 

As was discussed in Ref. [4JJ, we fix the coupling constant go* d-k from the partial width for the decay of D* — » Dtt 
which leads to go^Dix — 17.9. gwD'-n is related to go*D^ via HHChPT [56| . i.e., 

giro** = -: 9D * D * = 2-31 GeV" 1 . (19) 

In Ref. [57| . the HHChPT for the interaction of charmonium with open charm mesons was constructed. From this 
theory, one can yield the following relations for the coupling constants 



1 \/ m r>' 

gD*Dip n = ~7;gDD4> n T=, (20) 

2 Tn^yfrriD 

g X c DD = 3—^-g X coD-D* = -2v / 35i m r>v /m Xco i ( 21 ) 

g X dD*D = V2giy/m Xcl m D m D * , (22) 

g X c2D'D* = 2g 1 m D *^/m Xc2 , (23) 



with the expression for g\ as 



where f Xc0 is defined via the relation (0\cc\x c o(p)} — /x c o TO xco an d the QCD sum rules yield f Xc0 =510 MeV J58|- So 
that we have g 1 = -2.09 GeV" 1 / 2 . 

The electromagnetic interaction Lagrangian applied in our calculation can be decomposed as 

-Com = + £-kin gC + £-D*Di/>„7 + Co'D-y- (25) 

The effective Lagrangian £^ is yielded by substituting C9 V in Eq. © with Cfy s au s e with the explicit form 
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where the Wilson's line I(x, y, P) is defined as 



I(x,y;P) = / dz^(z). (26) 
Jy 

In the following calculation, the nonlocal vertex with one photon is necessary. This vertex comes from two sources: 
One is from the covariant derivative and the other one is from the expansion of the Wilson's line. One can easily derive 
the Feynman rule for the nonlocal vertex with one-photon coming from the covariant derivative. But to derive the 
Feynman rule for the photon from Wilson's line, one may use the path- independent prescription suggested in (59l - [6ll | . 
£^ ugc is from the gauged kinetic terms of the charged constituents D* + and D + , i.e., 



£g^ge = ie A fl (D-d fl D + ) + ieA fi [-D* ot -d^D* + a + -D* a -d a D* + ^ + -D*-^d a D* + a ]. (27) 



l -D* a -d a D* + »+ 1 -, 

And £D*Dip n -y is obtained by gauging the derivative coupling of the D*Dip n vertex in Eq. (|T5|) 

Cd-d^ = g D 'D^ a ^n:^D* a pD+ +H.c, (28) 

with £>*£ defined via the covariant derivative as -D*^ = (d + ieA) a D*p~ — [d + ieA) pD* a ~ . 
The effective Lagrangian Cd*Dj can be generally expressed in the form 

Cd- Di = eg D , +D+1 e^D+D;-F aP - eg D ,o D o^ ua ^D D;°F aP . (29) 

From the partial width for the decay D* — > one can determine the coupling constant gD*D-y- Using the data 
T(D* + — > D + j) — 1.54 KeV, we have g n * + n + 7 ~ 0.124 GeV" 1 . With respect to the isospin invariance, the decay 
D*° — y Z?°7r° can be related to D* + —> D°tt + [Tj|. And considering the fact that the branching ratios for the decays 
D*° — > Z?°7 and D* Q — > D°ir° are well measured [5(|, one can yield the coupling constant gD- a D- t — 0-51 GeV -1 . 
The relative sign between these two terms are fixed as in Ref. 62j. 



III. STRONG AND RADIATIVE DECAYS OF X(3872) 
A. Radiative X(3872) -> jj/ip and X(3872) -»■ -yip{2S) decays 



In this part, we study the radiative decays of A (3872). Taking into account the quantum numbers of the initial 
and final states, one can write the most general form of the gauge invariant matrix element for the electromagnetic 
decay of A(3872) -> j(qi)V(q 2 ) as 



iMxv-) 



ieMj^e^(p)e*(gi)e£(g 2 ) 
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(30) 



where e^ v is the polarization tensor for the tensor meson A(3872) satisfying the transverseness condition p^tpvip) = 
and the traceless condition e^(p) — 0. q is defined as q — q\—q2 and qi satisfies the transverseness condition qi-e(qi) = 0. 

In terms of the factors o-xy 1 defined in Eq. (|30[) . one can express the partial width for the radiative decay of 
A (3872) -> -fV as 



XV-f 



10m 2 , I ^ 



^XVi^XV^ XVi 



(31) 



where P* is the three-momentum of the decay products. The coefficients C^ v are functions of m x and my- We 
give their explicit forms in Appendix [SJ To get this partial width, we have applied the sum of the polarization vector 
for tensor meson 16311 



^ ] e HlVl{p) e u 9 v 9 .{p) — 2 (^1/^2^1 1^2 + ^tl^^l/^J g^Ml^l^, 



M2^2 J 



(32) 



polar 
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where 9^ u = -g^ + (p^p^/m^. 

In our present work we study the radiative decays of X — > jJ/tp and X — > jijj(2S). The relevant diagrams are 
illustrated in Fig. [21 After standard computation, we can yield the numerical results for the partial widths. Both of 



7 7 




D*-(D*+) D*~(D*+) 
(E) (F) 



FIG. 2: Diagrams contributing to the radiative decays of X (3872) — > 'yJ/ip, yt[>(2S). 
them depend on the size parameter A x and the mixing angle 4>, or equivalently, 9. 

B. Strong X(3872) -> J/iph decays 

Next, we study the strong decays of X(3872) — > J/tph with h as light pseudoscalar mesons tt + tt~ or ir + ir~ir°. 
Measurements of the invariant masses of pions indicate that these processes happen through X(3872) — > J/ipp for 
the X(3872) ->■ J/^tt+tt- decay ^MM but ^(3872) J/^uj for the X(3872) J/^tt+tt-tt decay So that 

in the following calculation, we only take into account the decays via light vector meson exchange. Explicitly, we 
consider the processes X(3872) — > J/ipp° — » J/ipir + ir~ and X(3872) — > J/tpw — > J/ipir + Tr~ir Q . Generally, one can 
write the matrix elements of these decays as 

iM{X^J/i>h) = iM a (X ^ J/W)- t^—- ^(9**-^ir L )M' 7 (V^h), (33) 

q(-my+im v T v m z v 

where V is the intermediate light vector meson with m v and Ty as its mass and total width respectively, qi is the 
momentum for the intermediate light vector meson V. 



G 



M a (V — > h) stands for the contribution from the V — h interaction which relates to the matrix element for V —> h 
decay via 



iM(V -> h) = ie a M a {V -> h), 



(34) 



with e CT as the polarization vector of the vector meson. 

The matrix element M a (X —> J/ipV) relates to the matrix element for the decay X —> J/tpV. Taking into account 
the transverseness of the polarization vectors and symmetric property and traceless of e^, we write the general form 
of M a (X -> J/tpV) as 

iM a {X^J/W) = iM^(m|,m^,« 2 )c^(p)eJ(ga) 



+G 



XipV 



(3) 

XipV 
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XipV 



* apfS Voq P q v q a + 0* f) 



}<v(p) e iste)i 



(35) 



where in this expression, p is the momentum of X(3872) and q = q\ — qi with q\ and q^ as the momentum of the 
intermediate light vector meson V and the final vector meson J / ip respectively. To write down this general form, we 
have considered that both M ti (p a — > tt + it~) and M^(uj — > 7r + 7r~7r°) are transverse to q p . In our calculation of the 

form factors G^ v (i = 1, 2, . . . , 7) we have taken into account the diagrams illustrated in Fig. [3] 
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FIG. 3: Diagrams contributing to the strong decays of X(3872) — > J/ipu — > J/ipit + -K ty°, J/ipp — >• J/ipir + ir . 

Using Eq. (|33|) one can write the differential decay rate into J/iph as a function of the invariant mass qi of the 
hadron h as 

where the physical region of q\ is rrih < 9i < Tlx — m-v anc ^ the coefficients Cij are given in Appendix IB II j dX\F(V — ¥ 
h) is the Lorentz invariant phase space integral for V — > h decay. For the p° — > 7r + 7r~ decay we have 

dUF (p° 7T+7T-) = -^-(fl? - 4m*) 3 / 2 , (37) 

Z47T (/I 
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with g p7rlT determined by the partial width for p° — > ir + n decay as 

g p ^ = 6.0 ± 0.03. 



(38) 



In the case of ui — > 7r + 7r~7r° decay, this process has been studied in Refs. [6J-|66( using the hidden local symmetry 
formulism. But in this formulism, the matrix element depends on two unknown parameter combinations, c\ — C2 — C3 
and C3, explicitly, 



iM»{uj -> ir + Tr-ir Q ) = -g 



16tt 2 F3 



e^q ;,q + ; a q-,p 



x {3( Cl - c 2 - c 3 ) + 2c 3 [D p ((q + + c?_) 2 ) + D p ((q_ + q ) 2 ) + D p ((q + <z+) 2 )] }, (39) 

where g = 5.80 ± 0.91 is the coupling constant introduced in the gauge of the hidden local symmetry and D p (q 2 ) is 
the propagator for p meson which is normalized as D P (Q) — 1 and at the leading order 



irr 



2 2 ' 
nip — q z 



(40) 



Since it is a formidable work to determine these two combinations in this paper we use the values determined in 
Ref. [17( which are provided in Appendix [Cj 



C. Strong X(3872) -> X cJ^° decays 



The decays X(3872) — > XcJ 71 " (J = 0, 1, 2) are important and the measurement of these processes can discriminate 
between the charmonium and the molecular interpretation of X(3872) in the case of J PC — 1 ++ [13] ■ We expect that 
this holds for J PC = 2 h case. In this subsection, we devote to study these processes. Generally, the matrix elements 
for these processes can be written as 



iM(X -> Xcott ) 
iM(X -> XdA 

lM(X Xc2TT°) 



Peq T e^(p)e* a (q), 



(41) 
(42) 



+G 



(3) 



(g^q v + g™q»)pP + (a o /?)] Je^e*^), 



(43) 



where in Eqs. (|42"T) and (j43|) q is the momentum of %ci and x C 2 but in Eq. (|4Tj) q = qi — q2 with qi and 52 as the 
momenta of the final states. In our explicit calculation of the effective coupling constants Gx XcJI », we include the 
diagrams illustrated in Fig. [4] 




FIG. 4: Diagrams contributing to the strong decays of X(3872) — > XcJi" • 
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With the general expressions for the matrix elements, one can express the general forms of the partial widths as 
T(X^ Xc0 n°) = |g^f A 2 K,< o ,^ )P x * c0 , (44) 

-i2 



T(X^ Xcl «°) = _^|A 2 K ! m^ 1 ,m 7r o)P x * ci , (45) 
1 3 

rpC^Xc97r°) = — =-[ V C% Q ]P* , (46) 

V ' 407Tm 2 L / ' ^Xc27T U XXc2^ U XXc2^' J Xc2' v. v 



»<j'=i 

where A(a 2 ,& 2 ,c 2 ) = (a 2 — (6 — c) 2 )(a 2 — (b + c) 2 ) is the Kallen function and P* cJ is the three momentum of the 

its rriir^n in A nT^n^i-^ IT3 01 



decay products. The coefficients C % 1 „ is given in Appendix IB 21 



IV. NUMERICAL RESULTS AND DISCUSSIONS 
A. Parameter fitting 

In the present analysis, because of the uncertainty of the V(3872) mass measurement, we take the binding energy 
defined in Eq. (fl"4"l) as AE = 0.5, 1.0 and 1.5 MeV. Then in our model, there are two parameters, i.e., the mixing angle 
(j) between the isospin singlet and isospin triplet components and the size parameter A x . In the numerical calculation, 
since we have no prior information from QCD, we take these two parameters as free ones to fit the following data 



B(X -> tK2S)) 
B(X -> jJ/ip) 

ff(V(3872) -> J/^TT+TT-TT ) 

B(X(3872) ^ J/iiw+n-) 
£(AT(3872) -> yj/i/)) 
B(X(3872) -> J/^tt+ti-) 



3.4 ±1.4 [|, (47) 

1.0 ± 0.4(stat.) ± 0.3(syst.) 0, (48) 
0.14 ±0.05 Belle 0, 

0.33 ±0.12 BaBar [8]. (49) 



We scan the mixing angle <f> from 0° to 90° with a step of 1° and the size parameter A x from 0.1 to 3.0 GeV with a step 
of 0.1 GeV. As an illustration, we plot the mixing angle <fi dependence of the branching ratio T(X —> r yip(2S))/T(X —> 
7J/V) at three typical values of A x = 2.4, 2.7 and 3.0 GeV in the case of AE = 0.5 MeV in Fig. [SJ This figure 
shows that any values for <j> cannot yield the ratio for A x = 2.4 GeV case, while the range for is constrained for 
A Y = 2.7, 3.0 GeV. In this way we determine the range of the two parameters using the data in Eqs. (|47I48I49I) . 

In Table. HI our results of these two parameters which can reproduce the data within la together with the values 
of relevant branching ratios are given. As can be seen, we can reproduce the data in the first two columns and the 
Babar data in the third column quite well but we cannot get the result of Belle. Also listed in Table. Q] are the values 
of the ratio B(V(3872) -> 7^(25))/B(X(3872) -> J/^n + n~) which was measured as 

B(X(3872) jip(2S)) „ „ n 

= 1.1 ±0.4, BaBar [|. (50) 



fi(V(3872) J/tJjtt+it-) 

Taking into account this ratio, one can see that the smaller values of the mixing angle <f> are preferred. So in the 
following, we restrict ourselves to the case with = 9°. 

In Table. [TT] we present our numerical results of the relevant coupling constants g x and g x . From these results 
we conclude that the effective coupling constants are stable against the size parameter A x . Once comparing with 
Ref. [32j . one may find that our results for the coupling constants are bigger than the corresponding ones given there. 
This can be interpreted as the following: In our case, A(3872) has quantum numbers J PC = 2 ^ so the coupling 
between AT(3872) and its constituents DD* is in P— wave. Compared with the case that V(3872) has quantum 
numbers J PC = 1 ++ , it needs stronger attractive interaction to compensate the repulsive interaction induced by 
angular momentum. 

Our result of the small mixing angle implies that the isospin singlet component is dominant. Explicitly, one can 
write the wave function of V(3872) in terms of the isospin basis as 

|X(3872)) = 0.988 x | A" (3872)) J=0 + 0.156 x |X(3872)) I=l . (51) 



9 



20 - 


1 ' 




' 1 


' 1 1 1 


- 


20 








A x = 2 .4 G eV 

A = 2.7GeV 

X 




15 J 








A = 3 .0 G eV 

X 




15 


10 - 




/ 

/ 








10 


5 - 


• I 


/' 


'.\ 

\\ 
' i 






5 




•' ,' 1 

.' /' '\ 


















- 














-I 1 1 1 1 1 1 1 1 1 

20 40 60 80 100 



M ikiig ang]s (j> 



FIG. 5: Mixing angle cj> dependence of the branching ratio V(X -> jtp(2S))/T(X -> jJ/ip) in the case of AE = 0.5 MeV. The 
grey band corresponds to the data with the error bar. 



TABLE I: Fitted parameters 9 and A x and the corresponding branching ratio. 



<t> Ax(GeV) 


A£(MeV) 


B(X^~fi>(2S)) 


B(X^ 






B(X^ 1 i,(2S)) 




B(X->J/i/.7r+jr~) 


B(X^J/i/.ir+ir-) 


B(X->J 7l/>7T + 7I— ) 


12° 2.6 


0.5 


2.089 


0.483 




0.218 


0.455 




1.0 


2.101 


0.461 




0.214 


0.450 




1.5 


2.113 


0.439 




0.210 


0.444 


11° 2.7 


0.5 


2.628 


0.591 




0.211 


0.555 




1.0 


2.640 


0.566 




0.208 


0.549 




1.5 


2.647 


0.542 




0.206 


0.545 


10° 2.7-2.8 


0.5 


2.158 - 3.007 


0.748 


- 0.740 


0.320 - 0.233 


0.690 - 0.700 




1.0 


2.165 - 3.014 


0.716 


- 0.705 


0.316 - 0.229 


0.684 - 0.690 




1.5 


2.168-3.017 


0.685 


- 0.673 


0.312 - 0.226 


0.676 - 0.682 




0.5 


2.430 - 3.486 


0.964 


- 1.035 


0.369 - 0.283 


0.897 - 0.987 


9° 2.8-3.0 


1.0 


2.433 - 3.504 


0.917 


- 0.971 


0.363 - 0.274 


0.883 - 0.960 




1.5 


2.433 - 3.516 


0.874 


-0.914 


0.358 - 0.266 


0.871 - 0.935 


data 




3.4 ± 1.4 [8] 


1.0 ±0.4 ±0.3 [7] 


0.14 ±0.05 [7J 


1.1 ±0.4 [8] 












0.33 ±0.12 [8] 





From this one might think that decay X — > J/ipir + ir~ would be strongly suppressed compared to X — » J/ipTr + ir~Tr° 
decay. However, this is not the case. One can understand this as the following: Concerning the fact that the mass of 
(3872) is around the threshold of J /ipV and following Ref. [l7[ one can simplify the matrix element (|35|) into the 
approximate form 

iM a (X J/i>V) = iG x ^v^ aap Paq v e^{p)e* 3 {q2)- (52) 

So that, in terms of the effective coupling constants Gx^u and Gxip P the ratio B(X — » J / ijm + ti~ k ) / B{X — > 
J/ipTT + w~) can be expressed as 

BjX -> J/^tt+tt-tt ) _ \G X ^\ 2 I d!l(X -> J/jiir+ir-ir ) 
B(X -> J/ifm+ir-) \GxtpA 2 J dH{X ^ J/^tt+tt-) ' 1 ' 

where J dH denotes the phase space integral. If one naively use the wave function (|5ip one can get \Gxxpui\ 2 l\Gxtpp\ 2 — 
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TABLE II: Effective coupling constants g x and g x from the fitted parameters <j> and A x . 



<j> A A .(GcV) 


A£(MeV) 




9° x 






0.5 


16.29 - 15.84 


12.64 - 


12.40 


9° 2.8-3.0 


1.0 


16.38 - 15.94 


12.68 - 


12.44 




1.5 


16.46 - 16.04 


12.72 - 


12.48 



TABLE III: Fitted parameters <j> and A x and the corresponding branching ratio related to the data. 







AE 




r(x^7vK2S)) 




IXJf ->■ J/^TT+TT") 




(GeV) 


(MeV) 


(KeV) 


(KeV) 


(KeV) 


(KeV) 






0.5 


2.085 - 1.872 


5.066 - 6.525 


5.447 - 6.842 


5.648 - 6.612 


9° 


2.8 - 3.0 


1.0 


2.082 - 1.864 


5.065 - 6.533 


5.253 - 6.606 


5.730 - 6.802 






1.5 


2.079 - 1.859 


5.058 - 6.537 


5.067 - 6.380 


5.801 - 6.977 



(0.988/0. 156) 2 ~ 40. But this ratio is strongly suppressed in case one includes the phase space factor, i.e., J dH(X — s- 
J/^tt+tt-tt )/ / dU(X ->• J/inr+n-) ~ 0.0127 which leads to B(X -> J/iPn+n~Tr°)/B(X -> J/^Tr+Tr") ~ 0.5. More 
precisely, without including the difference arising from the loop integrals in Fig. |3l one can yield 

~ . 9_ddu = (M) 
Gx^p 9 X - 5x 9dd p 

where we have used the numerical results of Table. ITT1 and Eq. (|16|). From this we yield B(X J j 'ipir + n~ ir°) / 'B(X —> 
J/'ipir + Tr~) ~ 1.0 which is consistent with the experimental data given in (J3HJ) • 



B. Partial widths for decays with J/ip or ip(2S) in the final states 

In Table. IIIII we give partial widths for decays with J/ip or i/j(2S) in the final states. From this we see the 
partial widths are of order of KeV. These partial widths have been computed in the case that AT (3872) has positive 
parity 0, H IH El . For the strong decays AT(3872) ->■ J/ipir + it- and AT(3872) J/^7r + 7r"7r°, it was found that 
the decay widths are around 50 KeV which are both about one order larger than our present results. But for radiative 
decays, the results depend on model closely. In the case that X(3872) with quantum numbers J PC = 1 ++ is a 
bound state of mesons 18], the radiative decay widths are of order KeV which are at the same order as our present 
results, while if X(3872) is a mixing state of molecule (without the charged DD* components) and charmonium 
component [Hj], the decay width of X(3872) — > jJ/tp is found to be of 100 — 200 KeV which is much larger than that 
computed from the hadronic molecule assumption. In Ref. [33| . an admixture model of molecule with the charged 
DD* components and the charmonium components was applied to study the radiative decays, the similar results as 
ours were yield. 

The inclusion of other components as was done in the case that AT (3872) has quantum numbers J PC = l ++ may 
change our results. One possibility is to include J/tpui and J/ipp in the wave function j26], [32J. This may increase 
the magnitudes of the strong decay widths and the results depend on the probability of J/ipui and J/ipP hi A^(3872) 
even in the case that only the long distance effect is considered. In this sense, if the strong decay widths for tensor 
A^(3872) are observed bigger than our present results, one may conclude that the tensor A^(3872) cannot be a pure 
DD* molecule and other component should be included. Another possibility is to regard A" (3872) as a mixing state 
of cc and DD*. One may borrow the lesson from the X(3872) with 1 ++ case [18| to naively expect that this change 
of the wave function of A'(3872) may improve the magnitude of radiative decay width of A"(3872). Of course, in case 
other constituent is included, the probabilities of the relevant components should be fitted from data again so the 
magnitudes of the partial widths should be studied in detail. 

From these discussions, we would like to stress that the precise measurement of the strong decay widths can provide 
some clues on the structure of A' (3872). 
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TABLE IV: Partial widths for X(3872) ->■ Xcjn decays. 






Ax 


AE 


r(x ->. xcotv ) 




T(X -> Xc27V°) 




(GeV) 


(MeV) 


(KeV) 


(KeV) 


(KeV) 






0.5 


22.41 - 21.97 


0.294 - 0.276 


207.5 - 8.335 


9° 


2.8-3.0 


1.0 


22.60 - 22.40 


0.296 - 0.281 


211.8 - 8.618 






1.5 


22.78 - 22.76 


0.296 - 0.285 


215.9 - 8.880 



C. Predictions of the X(3872) — > Xci^ decays 

As was discussed in subsection IIII CI we expect that the measurement of the decays X(3872) — > Xcj^ ( J = 0, 1, 2) 
can discriminate between the charmonium and the molecular interpretation of X (3872). With the expressions given in 
subsection IIII CI and the parameters fixed in subsection IIVAI the decay widths of X(3872) — > XcJ^ can be calculated 
explicitly. Table. ITVl is the summary of our numerical results. 

From this Table, one may find that the partial width of decay AT (3872) — > Xc2""° has a strong dependence on the 
size parameter A x . This can be understood as follows: From dimensional analysis, in Eqs. (I41I42I43I) . only 27r0 , 
which gives the dominant contribution to this process, has a positive mass dimension so it strongly depends on A x . 

From Table. IIVI one can yield the following ratio of the partial widths 

T(X -> Xc0 7r°) : T(X -> X ci^°) * 1:0.013, (55) 

which indicates that, compared with the decay X — > XcO 71 " , X — > Xci^° is strongly suppressed. From the expressions 
for the partial widths, using the expressions for the coupling constants in Eqs. (|19|21I22|) . one naively has the ratio 

T(X -> Xc0 7r°) 8 /4V2 \2 8 { 8V2 \* f 

r , v : - t-2"( —pr(m D . +3m D ~ -—[—=-m x ) ~ 110, (56 



which can approximately explain the ratio (|55[) . 

Similar to the discussions given in the previous subsection, the introduction of other components may change 
our numerical results. In case of the cc component is included, the magnitudes for the partial widths might be 
changed but the ratio for the partial widths must be kept since the cc component is a definitely isospin singlet so 
it does not contribute to the isospin violating decays AT (3872) — > Xc/tt - The same conclusion can be drown if 
J/tpui are constituents of X(3872). However complication arises if AT (3872) has a J /ipp component as was done for 
J PC — 1 ++ [Hj]. This is because the component J/^p gives a contribution to the isospin violating AT(3872) — > Xc.jk 
decays. 

From this, in the case AT(3872) with J PC = 2 the strong suppression of the decay AT(3872) — > Xci^° compared 
with the decay AT (3872) — > Xco^ ma y signal the pure DD* molecular structure of AT (3872). 



V. CONCLUSION 



In this paper, by regarding the hidden charm state ^(3872) as a DD* bound state, we studied its radiative and 
strong decays in the effective Lagrangian approach. We find, with an approximate probability of 97.6% isospin 
singlet component and a 2.4% isospin triplet component, the ratios B(X ji[)(2S))/B(X — > "fJ/ip) by Babar 0], 
B(X -+ J/ipTr + Tr-Tr°)/B(X -> J/^tt+tt-) by Belle and B(X -> -yJ/ip)/B(X -> J/^Tr+Tr") by Babar Q can 
be explained consistently, but the ratio B(X — > jJ/tp)/B(X — > J /%1)tt + it~) by Belle @ cannot be reproduced. We 
would like to stress that, in case of A"(3872) as a tensor meson, few percent of isospin one component of X(3872) 
in the wave function can accommodate the large isospin violating partial width since the phase space of the decay 
AT (3872) — > J/if)ir + ir~ir° is about one percent of that of the decay X(3872) —> J/^m + ir~ . 

With respect to the data, we fixed both of the parameters <j> and A x in our model, so that all the partial widths can 
be calculated explicitly. We calculated the strong decays X(3872) -> J/^tt+tt-tt , AT(3872) -> J/tfjir+n- , AT(3872) -s- 
Xcjtt {J = 0, 1, 2) and radiative decays X(3872) jJ/ip,-fip(2S). For the strong decays AT(3872) -> J/^tt+tt^ and 
AT (3872) — > J/ipTT + ir~Tr° , our results are several KeV which are both around one order smaller than the case A" (3872) 
with J PC = 1 ++ . For radiative decays, our results are both of order KeV. They are at the same order as those in the 
case that X(3872) with J PC = 1 ++ is a bound state of mesons [Hj], while much smaller than those obtained in the 
case that X (3872) is a mixing state of molecule and charmonium component [32l |. 
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For the X(3872) Xc/tt decays, we found that, compared with the decay X(3872) ->■ x c0 7r°, the X(3872) xa^ 
decay is strongly suppressed. The experimental observation of this suppression may be a signal of the pure DD* 
molecular structure of X(3872). 

Concerning the lessons from the study of X(3872) with J c = 1 ++ , we naively expect that the inclusion of 
other constituent in the wave function may change our numerical results of the relevant partial widths. Explicitly, the 
inclusion of J/tpu or/and J/tpP components may increase some partial widths for the strong decays while the inclusion 
of cc component may increase the partial width for radiative decay. As was noted, our above calculation cannot yield 
the ratio B(X —> / yJ/tp)/B(X —> J/ipir + ir~) by Belle [7] although we reproduced this ratio by Babar [8J. If the Belle 
value is preferred in the future experiment, in the molecular interpretation, we should include other constituents. We 
leave this inclusion of other constituents in our future work. 

Taking into the quantum numbers of X(3872) and its constituents, the coupling between X(3872) and DD* is 
via P— wave. We would like to mention that, although the meson exchanging model calculation indicates that it is 



difficult to form a DD* molecule via P— wave coupling [31|, we think that it deserves systematically and seriously 



analysis of the exchanged hadrons and coupling channel effects. In the present analysis, we do not consider origin of 
the binding force but regard X(3872) as a P— wave bound state of DD*. 

In conclusion, by regarding that the hidden charm meson X(3872) with J PC = 2 h is a DD* P— wave bound 
state, we found that our model with dominant isospin zero component can explain the existing date quite well. We 
also predicted the partial widths for X(3872) -> 7 J/^, X(3872) ->■ jip(2S), X(3872) -> J/^tt+tt-, X(3872) 
J / n~ tt° and X(3872) — > Xcjw ■ Comparison of these values with the future experiments will shed a light on the 
nature of X (3872). 

Appendix A: Explicit forms for the functions C"^^ for radiative decays. 

Here, we present the nonvanishing coefficients Cjc^^ for the radiative decay of X(3872) —> jip n with ipi = J/tp 
and ^2 = ip(2S): 

CxV, = 3w / m 2 M3A + iOmlmlJ, = -^A, 

C x% n7 = ,2, M3A + 10m 2 Y m^ - lOrn^J, C% nl = --^-^(m\ + m 2 JA 2 , 

r lb - 8 \ 2 r 22 -—x r 23 - 16 <m 2 m 2 )\ r 25 - 8 .\ 2 

°X-0„7 ~ 3m 2 A ' °X'*„7 ~ 3 ' ^V^7 - 3m 2 ^ m x °X^7 ~ 3m 2 A ' 

X x x 

C %m = o / 2 ( 3 ^x + 4m 2 JA 2 , C% nl = Ao~(m 2 x - m 2 JA 2 , 

r 35 _ 16 ( m 2 m 2 \>2 r 44 1 \3 ,-.55 _ 4 \3 / a 1 \ 

Cx ^ ~ 3< ( x " ^ jA ' Cx ^ ~ ' x ^ ~ 3< A ' ( > 

where A = (m 2 — m ) 2 . 

Appendix B: Explicit forms for the functions C lj for strong decays. 
1. Explicit forms for the functions C^j^y for X(3872) — 5* J/^fih decays. 

In this subsection, we list the nonvanishing coefficients C XJ ^ V for the strong decay of X(3872) — » J /iph: 

C xJ/^v = A(3A + 40m 2 m 2 M ), C% JH>y = ~X, 

C *W = 3TO 2 2 TO 2 A(3A + 20m 2 TO 2 M ), 

C xj/W = + 20m 2 - 2Om 4 J/0 ), 



13 



r<i5 

^XJ/ipV 



r <22 

^XJ/ipV 



^24 



n 26 



^xj/iiv 



^34 

^XJ/if/V 



r<35 

^XJ/i>V 



^f36 

^XJ/iiV 



f-,44 

^XJ/^V 



^xj/iiv 



^XJ/i>V 



\2i 2 2 i 2 \ 

2 ^2 A ( m x - m v + m J/v)> 



C 



17 



-A 



A(3A + 40m 2 m 2 ), 



A(3A + 20m 2 m 2 /v , - 20m 4 



r< 23 — ^ 



A(3A + 20m 2 m 2 ), 



\2/ 2,2 2 \ 



C 



27 



3m 2 to 2 m 2 T , , 

X V J/lp 



XJ ^ V - 3^ 
3(m 2 / + m 2 / ^)A + 40m 2 m 2 m 2 /j/ , , 



(m 2 -m 2 ^)A 3A + 40m 2 m 2 ^ 



m^ m j, , 



I 2 2 i 2 \ \ 2 

-(m - m +m JM )A 



C 



:S7 



1G 



XJ /^ " 3m 2 A ' 



1 



3m 2 (m 2 +to 2 , /v ,)A- 24m 2 m^A + 40m 2 m^m 2 - m 2 ,^) 2 



-(m 2 - m 2 - 3m 2 7/ JA 2 , 



,o46 



16 



(m 2 - m 2 7/j/) )A 2 , C X ji w - rn . 2 ni . 



\ 3 s~i66 
A ' ^XJ/ipV 



2 m 2 x 
1 



(m 2 - 3m 2 - m 2 /M )A 2 , 



A d , C 



77 



A 3 , (Bl) 



3m 4 v-v -j/*/". "«/^ m 2 m 2 "' "^/V" m 2 m 2 "' 3m 4 

X X J ftp XV X 

where A = \{m? x , m 2 ^, q 2 ) = (m 2 x — (fftj^ — <Zi) 2 ^ (^ m x ~ ( m J/iA + 9i) 2 ) i s the Kallen function with qi = \fq\ 



2. Explicit forms for the functions C lj for X(3872) — s> Xc27r decay. 



Here, we present the nonvanishing coefficients C Xx 27r0 for the strong decay of X(3872) — > Xc2^°- 

1 , A ^.n 1 



✓-.11 



r<i3 



^-<22 

L 'XXc2 7I- 



^-<23 



,-r33 



576mt 



-A 4 , 



X""Xc2 



ri\2 



36mt-mt 

A Xc2 



A (A + 6m x mi. ), 



1 



36 "4 m Xc2 



(m^ 



X c2 _ TO ^o)A 3 , 

A 2 + 30m^m 2 A + 180m^m 4 



9 "4 TO Xc2 



9m 4 m 4 ^ + m ' c2 ~ m ' o)A(A + 15m * m L)' 



I -X'"Xc2 

1 



A 2 (2A + 17m 2 f m 2 ), 



where A = A(m 2 ¥ , m 2 c2 , m 2 ) = (m 2 x - (m Xc2 - m w o) 2 ) (m 2 x - (m Xc2 4-m^n) 2 ^ 



is the Kallen function. 



(B2) 



Appendix C: Decay widths for lo — > ty + -k ■k . 

From the explicit study of Ref. [13] we have the following expression for the cj — > 7r + 7r~7r° decay that is used in this 
work. 



dIlF(uj -J- 77+77^77°) 



(cos 9 V + \/2 sin V ) 2 
^Ffqf 



2 Aqi-m^Y 



im'i 



ds 



X>, 



^.2 12 2\2 

m 7r{ m 7T ~ Ql) -S12S23S13 
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^9 p7T7rC VV7r 



(l-i[/p(*12) + />38) + /p(»3l)]) 



where sx2, S23 and S13 are the invariant mass squares of the three pion pairs and 



, 2;5 



4m^ 
S12 



and 



s - M$ + im v Tv '' 



(CI) 



(C2) 



(C3) 



is a vector meson resonance factor that vanishes at s — 0. The pion decay constant is F„ = 92.4 MeV, the values of 
the parameters C v ^ and q C mm F 2 /m 2 , are given by Eqs. (|C4|) and (|C5[) . and the value of the light vector meson 
mixing angle 9 V is given by Eq. (IC6[) [17j 

(C4) 
(C5) 



+ 2 g(mw C vvw F*/m 2 v = (8.03±0.48)/(16^ 2 ), 
g p ^C m Fl/ml = (5.1±0.65)/(16^ 2 ). 

The vector meson mixing angle is given by (68l . l69j 

cos6>„ =0.51 ±0.01. 



(C6) 



The function cos0„ + \/2sin#„ used in the numerical function is cos#„ + \/2 sin V ~ 1.73 i 0.01. The errors in the 
parameters in Eqs. (|38p . (|C4[) . (|C5p . and (IC6|) arc determined using the uncertainties in the measurements of the 
vector meson decay widths only. 
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